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O ; Abstract 



We study the phase diagram of the one dimensional (ID) U — V model at 
' quarter filling in the most general case where the on-site and first-neighbour 

interactions U and V can be both attractive and repulsive. The results have 
been obtained using exact diagonalization of small clusters and variational 
techniques, as well as exact results in various limits. We have analyzed 



^ ■ four properties of the groundstate: i) whether it is insulating or metallic; 

ii) whether it is homogenous or phase separated; iii) whether it has a spin 
gap; iv) whether it has dominant superconducting fluctuations. With eight 
phases, the resulting phase diagram is unexpectedly rich. The four phases 
not found in the weak coupling limit are: i) an insulating phase when U and 
V are large enough; ii) a region of phase separation when V is attractive; iii) 
another region of phase separation when V is large enough and U small; iv) 
a region with dominant superconducting fluctuations when V is intermediate 
and U small. The actual nature of this last phase, which has pairs but no 
spin gap, is not fully clear yet. 
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I. INTRODUCTION 



The phase diagram of one- dimensional models of correlated fermions is well understood 
in the limit of weak interactions. In that case, the model can be mapped using perturba- 
tion theory and a cut-off procedure onto the so-called Fermi gas model (also referred to as 



interactions, perturbation theory cannot be used to perform such a mapping, and the deter- 
mination of the phase diagram is in most cases a very difficult problem. The exceptions are 
models that have an exact solution, like for instance the ID Hubbard modeli. In the case 
of that particular model, no qualitatively new physics appears in going from weak to strong 
coupling. So it is legitimate to wonder whether there is a real need of studying the strong 
coupling regime of more complicated models. We believe this is the case for at least two 
reasons. First, real materials whose properties are of one-dimensional character are not nec- 
essarily in the weak-coupling regime, and understanding their electronic properties requires 
some knowledge of the properties of intermediate or strong coupling models. A very impor- 
tant example is the case of the organic conductors of the (TMTSF)2X family. In a large 
temperature range, they exhibit low-energy properties that have a clear one-dimensional 
character and that are compatible with the Luttinger liquid picturei^ if the exponent a 
defined by N{u) — N{uf) oc Q{ujf — uj){uj — cu/)" is slightly larger than 1, where N^u) is the 
density of states. Such an exponent is impossible to get with weak interactions which always 
yield a <^ 1. The only hope to reach such a value is to turn to strong coupling models. Let 
us note that the Hubbard model is not able either to reproduce such an exponent: a is at 
most equal to 1/8, a value attained in the very strong coupling limit U = +0x3!'!. 

Second, it has been known for a long time that the strong coupling regime of extensions 
of the Hubbard model can give rise to new physics. For example, it was first shown by 



model yields a Metal-Insulator transition for V = 2t when U is infinite in the case of 
quarter-filling, a phenomena not observed in the weak coupling regime of the same model. 



(7-ology), a model whose properties are 



well establishedlllS. In the opposite limit of strong 





that adding a repulsion term V between nearest neighbours to the Hubbard 
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More generally, commensurability between the number of particles and the number of sites 
can give rise to a Metal-Insulator transition in the weak-coupling limit only at half-filling, 
while this effect can in principle occur for any commensurate filling if strong interactions 
are considered. 

In this paper, we concentrate on the properties of the U — V model at quarter filling 
defined by the Hamiltonian 



where U and V are allowed to take positive or negative values, i.e. to describe repulsive or 
attractive interactions. An account of the properties of this model in the case of repulsive 
interactions has already been givenil. Unexpected results showed up in the limit U small, 
V large, and an extension of the model to negative couplings was clearly needed in order to 
clarify the properties of the model in the large V limit. 

Before we start the description of the phase diagram of the model of Equation (|I|), let 
us recall what is meant by phase diagram in the context of one-dimensional system Jili. 
In general, phases are distinguished by some kind of long-range order. If one sticks to 
that definition, there are only three possibilities for ID systems: i) Insulating phase if 
there is a gap in the charge sector; ii) Metallic phase if there is no gap in the charge 
sector; iii) Phase separation if the ground-state is not homogenous. Any kind of magnetic 
or superconducting long-range order is forbidden in ID because of quantum fluctuations. 
However, the correlation functions do not decay exponentially but as power laws, that is very 
slowly, and it is useful to distinguish regions in the metallic phase according to the correlation 
function that decays most slowly. This is useful because any infinitesimal coupling between 
the chains, that will be present in real materials, will introduce long-range order associated 
with that correlation function. According to g-ology, there are at least four different types 
of metallic phases that can be distinguished by the value of the critical exponent Kp (to 
be defined below), and by the presence or absence of a gap in the spin sector. Quite 
generally, superconducting fluctuations dominate if Kp > 1, while charge and spin density 




(1) 
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fluctuations dominate if Kp < 1. Now, it tliere is no spin gap, triplet and singlet pairing 
correlation functions on one hand, charge density and spin density correlation functions on 
the other hand, have the same power law behaviour, and one has to look at logarithmic 
corrections to see that triplet pairing and spin density fluctuations dominate over singlet 
and charge density fluctuations respectively. If there is a spin gap, triplet pairing and spin 
density fluctuations decay exponentially, and singlet pairing and charge density will be the 
dominant fluctuations. To summarize, one can expect to find six types of phases: 

1) Metallic with dominant spin density fluctuations if Ac = = and Kp < 1; 

2) Metallic with dominant triplet pairing fluctuations if Ac = A^, = and Kp > 1; 

3) Metallic with dominant charge density fluctuations if Ac = 0, A^ > and Kp < 1; 

4) Metallic with dominant singlet pairing fluctuations if Ac = 0, Aj, > and Kp > 1; 

5) Insulating if Ac > 0; 

6) Phase separation if the compressibility of the homogenous phase is negative. 

The paper is organized as follows. In Section II, we summarize the results one can get 
from g-ologj in the weak coupling regime. In Section III, we present the general phase 
diagram and provide detailed explanations on how the various boundaries have been deter- 
mined. In Section IV, we give some exact results that can be obtained in the V = +oo 
case and that are useful to check the numerical determination of the boundaries. Finally, a 
detailed discussion of the large V, small U case is given in Section V. This includes a deriva- 
tion of an effective Hamiltonian in that limit and some numerical results of flux quantization 
used to analyze the presence of pairs in the ground-state. 

II. THE WEAK-COUPLING PHASE DIAGRAM 

When the interaction terms are small compared to the Fermi velocity, one can rewrite 
the Hamiltonian in momentum representation and map it onto the g-ology model. For the 
Hamiltonian of the Eq. (||) at quarter filling the g parameters are given by 

gi = U + 2Vcos2kFa = U 
5 
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U + 2V 







U + 2V 



94\\ 



2V 



(2) 



The coupling is equal to zero because we are away from half filling. 

The schematic phase diagram for the quarter-filled model is presented in Fig. |l|, where 
the different phases are numbered from (1) to (4). If f/ > ((^i > 0), we scale to the 
Tomonaga-Luttinger model. From the scaling invariant 2g2 — gi the fixed-point value of the 
forward scattering is g2 = U/2 + 2V. For U > —V/4 there are strong charge- or spin-density 
wave fluctuations (region numbered by (1) in Fig. P and for U < —V/A {g^ < 0) large 
superconducting fluctuations are present (phase (2)). On the other hand, if U < {gi < 0) , 
we scale to strong coupling. We know from the solution along the Luther-Emery line that 
there is a gap in the spin spectrum, and charge- density (3) or superconducting fluctuations 
(4) are favoured when U -|- V"/4 is positive or negative respectively. 



To determine the general phase diagram, we must calculate the charge gap Ac, the spin 
gap A^, the critical exponent Kp and the compressibility k. 



Let us start with the charge gap. It is given by A^ = lim^^^+oo A(L; L/2), where A(L; N) 
is defined by 



In that expression, Eq{L; M) is the ground-state energy of M particles on L sites and can 
be obtained by exact diagonalization of small clusters using Lanczos algorithm. Such a 



III. THE GENERAL PHASE DIAGRAM 



A. The charge gap 



A(L; N) = Eo{L- N + 1) + E^{L- N - I) - 2E^{L- N) 



(3) 
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procedure has been used by Spronken et a for the case of spinless fermions with nearest 
neighbour repulsion, and the critical value Vc = 2t could be reproduced quite accurately by 
fitting A(L; A^) with polynomials of 1/L. In our case, the Hilbert space is much larger, and 
we have results for 3 sizes only (L=8,12,16). So, we cannot expect a very good accuracy. In 
fact, we could not get meaningful results by fitting these 3 values with A + B / L + C/L"^. 
However, fitting any pair of values with A + B/L gives reasonable results that do not depend 
too much on the pair chosen to do the fit. Besides, as far as the large U case is concerned, 
the best result is obtained by fitting the results for L = 8 and 16. This is probably due 
to the fact that the 12 site system is quite different from the other two: To get smooth 
results as a function of L, one must choose the boundary conditions so that kp is one of 
the allowed k values, that is periodic ones for 8 and 16 sites, and antiperiodic ones for 12 
sites0. It turns out that the insulating region is confined to the U > 0,V > sector. The 
resulting boundary goes from {U,V) = {+oo,2t) to {U,V) = (4t, +oo). The first point is 
just the exact result obtained by Ovchinnikov some years agoi, while the second one is an 
exact result that we derive in Section III. 



B. The spin gap 

We now turn to the spin gap. It corresponds to the energy needed to make a triplet 
excitation from the singlet ground-state and can be calculated as = lim^^+oo As(L; L/2), 
where As{L; N) is defined by 

A,(L; N) = EoiL; N; = 1) - Eo{L- N; = 0) (4) 

Eo{L; N; = s) is the ground state energy of a system of particles on L sites in the 
sector 5*^ = s. As for the charge gap, and presumably for the same reasons, the best way of 
extrapolating the finite-size results is to do a linear fit of the results obtained for L = 8 and 
L = 16. The resulting boundary goes through the point {U, V) = (0, 0) with a slope parallel 
to the V it should according to the weak coupling results. Besides, our numerical 
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results seem to indicate that it ends at the point (f/, V) = (— 4t, +00), a result we also prove 
in Section III. 



C. The critical exponents 

The nature of the ground-state and of the low-lying excitations when there is no charge 
gap is a very difficult problem. In principle, one can get the excitation spectrum for finite 
systems and see what the elementary excitations are. Hopelessly, with the sizes that are 
available, it is not possible to distinguish any structure that will remain in the thermody- 
namic limit with any degree of confidence. The only thing that is possible at the moment is 
to assume that the spectrum has a certain structure and to check if this is compatible with 
the finite-size results. In the case of ID systems, we know that, in the weak coupling regime, 
the metallic phase is a Luttinger liquid. Noting that this is also true for the Hubbard model 
for all U, it is a reasonable hypothesis to assume that this is still the case for the model of 
equation (|l]). Then one can deduce the parameters of the model from the low energy part 
of the spectrumlii'lii obtained numerically for small clusters, and check a posteriori if these 
parameters are consistent with the hypothesis. One can estimate the ratio Up/ Kp by using 
its relation to the compressibility k 

TT^ ^ (r) 

2 Kp n^K ^ ' 



L / Eo{L; N + 2) + E,{L- iV - 2) - 2E^{L- N) 
nA 4 



(6) 



Equation is the finite-size approximation to the compressibility, Eq{L; N) being the 
ground-state energy calculated with suitable boundary conditions. The charge velocity can 
be obtained from the low-energy spectrum as 

Up = {E,m L/2- 5 = 0)- E,{L- L/2))/(27r/L) (7) 

where Ei{L; L/2; S = 0) is the energy of the first singlet excited state. The analysis of the 
t — J model by Ogata et allll was based on these equations. 
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In our study, we have also used another relation that holds for Luttinger liquid, namely 

ao = 2upKp (8) 

where ctq is the weight of the Drude peak of the conductivity. In ID systems, uo can be 
obtained simply from&i 



7r92£;o(0) 



L 



(9) 
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where Eo{(j)) is the ground-state energy as a function of a phase due to a flux $ = L0 
through the ring. Equations provide us with 3 independent conditions on Up and Kp. 

This is very useful for several reasons. First, the conductivity is much simpler to evaluate 
numerically than the compressibility, which involves systems with + 2 particles. Besides, 
we have good reasons to believe that equation @ does not give a rehable estimate of the 
compressibility when V is large, and it is important to be able to determine Kp without its 
help. 

It turns out that there are two disconnected curves Kp = 1. One curve goes through the 
point {U, V) = (0, 0) and corresponds to the boundary found in the weak-coupling case. The 
slope found at (f/, V) = (0, 0) in our numerical determination of this boundary is in good 
agreement with the one predicted by the weak-coupling analysis. There is however another 
line Kp = 1 in the small U, large V region which has no counterpart in the weak coupling 
case. This line goes through the point (f/, V) = (0, ~ 8t) and seems to extend toward the 
points (4t, +oo) and (— 4t, +oo). It is not possible to decide whether these points really 
correspond to the limits of this line on the basis of our numerical results however: The 
estimate of Kp becomes unreliable when V gets very large. As we couldn't get exact results 
about Kp in the V = +oo case, this behaviour remains a conjecture. Our numerical results 
also suggest that this line lies entirely in the region where there is no spin gap, although we 
cannot be totally sure either. 

To complete the analysis of the metallic phase, we now have to check the validity of the 
hypothesis we made about the excitations of the system, namely that they are those of a 
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Luttinger liquid. We can check this hypothesis in two independent ways. First, we can use 
the consistency of the three equations used to determine Kp and Up to check the assumption 
that the system is a Luttinger hquid: Equations (^|]) and (|^||) are typical of Luttinger 
liquids and will be violated if we have an instability to another phase (e.g. CDW insulator). 
A convenient way to measure the consistency of the three conditions is to calculate the 
ratio aQ/ixn^nu^p which equals 1 for a Luttinger liquid. We have performed a systematic 
evaluation of this ratio along lines in the {U, V) plane. In the region where equation can 
be used, i.e. when V is not too large, we found that, to a good accuracy, this ratio equals 1 
in the metallic phase, that it drops rapidly when one enters the insulating phase, and that 
the transition point was in good agreement with our previous determination of the phase 
boundary. 

Another way to check that we have a Luttinger liquid is to extract the central charge 
from the finite-size scaling of the ground-state energy densityEl 

eo(L) = eo(+oo) - ^^^^|^c + o(l/L^) (10) 

Ua was obtained from 

u„ = {Ei{L; L/2; = 1) - E^{L- L/2))/(27r/L) (11) 

where Ei{L; L/2; = 1) is the energy of the first excited state in the sector = 1. 
Comparing our results for L=8, 12 and 16, we found that this scaling form was accurate 
except when V is too large, and that the central charge c determined in this way equals 1 
within 2%, in good agreement with the exact value c = 1 for Luttinger liquids. 

To summarize, we have been able to prove numerically that, when it is metallic, the 
system is a Luttinger liquid, except when V is large, in which case finite-size effects prevent 
us from checking this hypothesis. 
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D. The compressibility 



Finally, we must calculate the compressibility to determine whether the homogenous 
phases considered so far are stable against phase separation. Reliable estimates of the 
compressibility are very difficult to extract in the interesting regions where it changes signs 
due to huge finite size effects. So we had to turn to variational techniques to estimate 
this quantity. More precisely, we have used an extended Gutzwiller ansatz to calculate the 
ground-state energy as a function of density. Details can be found in Appendix A. As for Kp, 
we found two disconnected lines of phase separation. One is entirely located in the < 
half plane and goes from ([/, V) = (+00, —t) to (f/, V) = (—00, 0). The first point is exactly 
the value obtained by other authors from the mapping of that model onto the anisotropic 



XXZ Heisenberg model in the limit U = +00 (see e.g. Ref. 0). The other limit is not 
surprising either: For U = —00, the groundstate is a collection of static pairs. If > 0, 
they will tend to stay apart, while if < 0, they will gain energy by staying as close as 
possible. The other phase separation line lies in the small U, large V region. It goes through 
the point (f/, V) = (0, 15t) and apparently extends toward the points (U, V) = (4t, +00) and 
(f/, V) = (— 4t, +00). In fact, we will prove in the next Section that these points lie on the 
boundary. So this variational approach to the determination of phase separation looks quite 
reliable. 



E. The phase diagram 

The resulting phase diagram is shown in Fig. |^. It contains eight different phases. The 
physical origins of most of them is already clear. The insulating phase is a spin density 
wave that occurs because of commensurability and was discussed in some details before!. 
The nature of the ground-state in the V < phase separation region is also quite easy to 
guess. If U is large enough, one will have a region of density 1 with one particle per site 
and an empty region, while if U is small or negative, one will have a region of density 2 
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and an empty region. The limit between the two cases is presumably the continuation of 
the As = line, although we have not tried to determine it. Besides, four of the metallic 
phases have a weak-coupling counterpart, and they simply correspond to cases 1) to 4) of 
the Introduction. The two phases that occur when V is large and U is small are not so easy 
to account for however. The rest of the paper is devoted to a detailed discussion of this 
region. In particular, the origin of the phase separation will be clear from the solution of 
the V = +00 case. 

IV. THE V = +oo CASE 

In this limit a pair of electrons on the same site cannot be broken, since it would mean 
hopping to the nearest site at an energy cost of V = +oo. The treatment of the model in 
this limit is greatly simplified by this fact. The unpaired fermions can move in the regions 
located between pairs. Moreover, unpaired fermions cannot be on neighbouring sites, so the 
spin is not relevant and its contribution is only to increase the degeneracy of the energy 
levels. In Appendix B an effective Hamiltonian is derived for the case where the V is much 
larger than U and t, and the V = +oo case is described by the Hq (Eq. pSj) . 

The eigenstates can be classified according to the number of pairs M, to their position 
im L, 1 < m < M) and to the number of unpaired electrons between pairs m and 
m + 1. The number of sites between pairs m and m + 1 is given by Lm = im+i —im — ^- Now 
the essential point is that the subsystem made of Nm unpaired fermions located between 
pairs m and m + 1 has the same energy levels as a system of Nm free spinless fermions on 
a finite chain of length — — 1. This can be shown in two steps: i) there is a one to 
one correspondence between the states, ii) the matrix elements of the Hamiltonian are the 
same after mapping. 

The reduction of size is due to the constraint imposed hy V = +oo, which implies that 
two fermions cannot be on nearest neighbour sites and thus that the number of available 
configurations is greatly reduced. Now, if we consider the entity made of a fermion and the 
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site right to it as a new particle, the number of available sites is reduced by the number of 
fermions between the two pairs and by 1 because of the pair at the boundary. For example, 
if we consider a state where there are 7 sites and two fermions between two pairs, there are 
6 configurations and the correspondence reads 

|cie I e J, eeede) \dsseed)' 
\de] ee [ eede) — > Idsesed)' 
\de t eee I ede) \dseesd)' 
\dee t e J, eede) — > \dessed)' 
\dee t ee J, ede) —>■ \desesdy 
\deee t e J, ede) — > {deessd)' 

where e, d and s stand for empty, doubly occupied and single occupied sites. Here the 
number of sites available for the free fermions is 4 instead of 7. 

Then it is easy to see that these new particles obey the Pauli principle, and that the 
action of Hq is to let these new particles hop with an amplitude —t whenever there is a free 
site next to it. For instance, 

Ho\de t ee J, eede) — —t {\dee t e J, eede) + \de ^ e I eeede) + \de t eee I ede)) (12) 

can be translated as 

H^ldsesedy = —t{\dessed)' + \dsseed)' + \dseesd)') (13) 

So the energy levels are those of Nm spinless fermions on Lm — Nm — 1 sites. 

Let us remember that if we have N spinless fermions on L sites with open boundary 
conditions, then the states can be enumerated by quantum numbers q — 1 . . . N with the 
energies — 2t cos qn/ (L + 1), so that the ground state energy is given by 

E,,{L,N) = -2tJ2cosy^ (14) 

( 7V + 1/2 \ / 1 1 

= t - d sin —t: sin tt (15) 

V L + 1 y V 2L+1 ) ^ ' 
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Thus, the lowest energy of the configuration above is given as 

M 

E^MU + Y^ E,i{Li -Ni-1, TV.) (16) 

i 

Now we need to get the configuration of doubly occupied sites for which the energy 
derived above is the lowest. One can easily check that if we divide a system of spinless 
fermions into subsystems, so that the numbers of sites and particles are conserved, the 
energy will always be higher than that of the original system. Note that this result holds 
because we are considering open boundary conditions. As a consequence the lowest energy 
is reached if all the pairs sit together. Hence, the lowest energy for a system having M pairs 
is given by 

E = MU - Ea{L -N,N - 2M) (17) 

where we have taken into account the fact that the effective space for the N — 2M unpaired 
fermions is reduced from L - 2M + 1 to L - 2M + 1 - {N - 2M) - l^L-N. 
In the thermodynamic limit (L — > oo) this is given by 

^/L = !^t/--(l-n)sin7r-^(l-m) (18) 
2 TT 1 — n 

where m = 2M/ N is the proportion of fermions forming pairs and n = N/ L is the density 
of the fermions. 

This energy should be minimized as a function of m. II U > Uc — — 4icos7rn/(l — n) 
the energy is minimized for m = and there are no pairs in the system. If —4t < U < Uc, 
the minimum is reached for m given by 

2n — 1 1 — n 1 / , X /.^v 

m = \ arccos(-C//4i) (19) 

n n TT 

Finally, when U < —At the energy is minimal for m — 1, which means that all the fermions 
are paired. 

For n = 1/2 (quarter filled system) Uc = 4:t. Above 4t every second site is occupied by 
one fermion with spin up or down. This would correspond to a filled band in our effective 

14 



model of spinless fermions. However, if we create a pair, it would cost us an energy [/, and 
we can gain —4t in kinetic energy, since we have created two holes at the top of the filled 
band each having energy —2t (see Fig. Thus the creation of pairs becomes favorable 
when U < At. This estimate is in complete agreement with the value of Uc obtained from 
Eq. I for the charge gap, where Eo{L,N) = (Fig. |a), Eo{L,N + 1) = U (Fig. |b) and 
Eo{L, N — 1) = —4t (Fig. |5c,d) and Ac = U — At becomes zero when U = At. This proves 
that (f/, V) = {At, +oo) is on the boundary of the metal-insulator transition. 

Let us now look at the spin gap. When U < —At we have only local spins in the ground 
state, and A^, > . When U becomes larger than —At, paired fermions appear in the ground 
state. But we saw above that the energy is independent of the spin, which means that 
As = 0. So {U, V) = {—At, +oo) lies on the boundary of the spin gap. 

Finally, it is clear from Eq. (|1^ that there is phase separation: The ground state consists 
of a domain of density 1 where all the pairs are packed and a domain of density < 1/2 where 
the unpaired particles can move. So the points {U, V) = {—At, +oo) and {At, +oo) lie on the 
boundary of the phase separation. 

V. DISCUSSION OF THE LARGE V REGION 

The nature of phase (8) is quite puzzling. On one hand, our numerical results indicate 
that this phase has no charge gap, no spin gap, and that Kp is larger than 1. Then, according 
to g-ology, the triplet pairing fluctuations should dominate. On the other hand, if we try to 
understand why this phase might be superconducting, our results concerning the V = +oo 
suggest that it is due to the presence of local pairs in the ground-state. The argument goes 
as follows: When V = +oo and —At < U < At, the ground-state consists of local pairs that 
sit together, and of unpaired particles that can move freely (apart from the V interaction) 
in the part of the sample not occupied by pairs. Then, when V decreases, the pairs start 
moving around. This can be seen by deriving an effective Hamiltonian in the limit of large 
but finite V (see Appendix B). The system clearly consists of two types of carriers, local 
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pairs and unpaired particles, that can hop, cross, or exchange. We have not attempted yet 
a careful study of this effective Hamiltonian. One scenario as to what happens when V 
decreases is the following. As long as V is not too small, the pressure exerted by the free 
carriers is sufficient to force the pairs to stay together, and there is still phase separation. 
However, when V is small enough, the kinetic energy gained by the pairs is comparable to 
that of the free carriers, and they can move freely in the sample. The ground-state is then 
homogenous and metallic. Superconductivity would then arise from the presence of locals 
pairs. But these pairs are located on a single site and must be singlet. Then this picture 
would favour singlet superconductivity, in contradiction to the prediction of gf-ology. 

Let us note that there is no inconsistency so far. On one hand, we have not been able 
to check that the system is a Luttinger liquid when V is large, so that the predictions of 
gi-ology shouldn't be taken too seriously in that region. On the other hand, we have no 
precise indication that there are still pairs in the ground-state when V is so small that we 
are outside the region of phase separation. In fact, an alternative scenario to that outlined 
above for the large V case is that the pairs disappear when one leaves the region of phase 
separation. Instead of having to reconcile these points of view, it looks more like we have to 
choose between one of the following possibilities: i) the system is a Luttinger liquid, there 
is no pair in the ground-state, and the triplet pairing fluctuations dominate; ii) the system 
has two types of carriers, local pairs and unpaired particles, it is not a Luttinger liquid, and 
the singlet pairing fluctuations dominate. How can we choose? 

The most natural thing to do would be to look at the low-lying excitations of the system 
to see whether there are of the Luttinger-liquid type or not. We have indeed tried that, but 
with no success. For systems as small as 16 sites, it is impossible to say anything by looking 
at the spectrum directly. There is however an indirect way of testing whether there are pairs 
in the ground-state by studying the response of its ground-state energy to a magnetic flux. 

To be more specific, suppose that a closed ring is threaded by a magnetic flux The 
units are such that $ = 27r corresponds to one flux quantum hc/e. This can be taken into 
account by replacing the kinetic energy by 
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-tY.i^^'M.c.^i. + h.c.) (20) 

where L is the number of sites. Then, for a noninteracting system the function /($) defined 
by /($) = Mvul^oo L{Eq{(^) — i?o(0)), where -E'o('^') is the ground state energy for a given 
flux, is periodic with period 27r. However, if the ground state consists of pairs, hke for the 
negative U Hubbard model, the function /($) is periodic with period vr. If the ground state 
has both types of carriers, then we expect to be in the intermediate situation where this 
function is periodic with period 27r but has a local minimum at vr. 

For finite systems, one finds almost systematically two minima. This is even true in the 
case of the Hubbard model with repulsive U . The only meaningful information is contained 
in the limiting behaviour when the system-size becomes infinite. With systems consisting of 
8, 12 and 16 sites, all we can see is the trend. Some results are given in Fig. ^. Quite clearly, 
when U is negative, but still inside phase (8), the local minimum at vr gets deeper when 
L increases, which suggests that there are pairs in the ground-state in the thermodynamic 
limit. For larger values of U , there is still a well defined local minimum, but it tends to 
decrease with the size of the system, so it is impossible to know whether it will survive in 
the thermodynamic limit. 

To summarize what we have learned about phase (8), it seems to us likely that when 
f/ < there are pairs in the ground-state, that the system is not a Luttinger liquid but 
has a more complicated spectrum including excitations of both the pairs and the unpaired 
particles, and that it has strong singlet pairing fluctuations. When U becomes positive, the 
nature of this part of the phase diagram is still very much an open question. There might 
still be pairs around, so that the whole region has the same behaviour, but the pairs might 
have disappeared as well, in which case the system would be a Luttinger liquid with triplet 
pairing excitations. More work is needed to settle this issue. 
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APPENDIX A: CLUSTER GUTZWILLER APPROXIMATION 

The ground-state energy has been estimated with the help of a variational wave function, 
namely a Gutzwiller-like wave functionll extended by intersite correlations. That intersite 
correlations might be essential was first pointed out by Kaplan alM, who showed that the 
correct —t^/U energy can be obtained in the case of the half-filled Hubbard model by includ- 
ing the intersite empty-doubly occupied correlation in addition to the on-site correlation. 
This kind of variational wave function has been applied to different models, and it proved 
to be useful to get results in regions where other methods have difficulties. 

The variational wave function we used can be written 



\^eg) = P\FS) 



(Al) 



where the operator P contains the intersite correlations and is given by 



The projectors A j at site i are defined by 



(1 - n|,i)(l - if/i = e; 

(1 -n|,i)n|,i, if/i=i; 

^T,i^i,i' if = 



(A2) 



(A3) 



The Ally's are the weighting amphtudes of the different nearest neighbour intersite corre- 
lations. For instance, if an empty site and a doubly occupied site are nearest neighbours, 
the weight of the configuration is multiplied by Xed- Altogether there are 16 different X^u^s, 
but they are not independent. Different symmetries, like A^i^ = A^/^ and zero magnetization, 
reduce the number of independent weighting amplitudes to 7. Furthermore, the density of 
particles being given, there are constraints to be satisfied in the minimization process. 

The difficult step is to calculate {'^eg\H\'^eg)- In general, this cannot be done 
analytically^, and the only "exact" way is to do it with variational Monte CarloB. There 
is however an approximate way of calculating this matrix element known as the cluster 
Gutzwiller Approximation!^ which leads to relatively modest numerical tasks, and which 
has been shown to give very good results on various models. In this scheme, a cluster of 
four consecutive lattice sites is chosen, and they are labeled by 1,2,3 and 4. On this small 
cluster we can have 4^ = 256 different configurations. The weight of each configuration is its 
weight in the Fermi sea multiplied by the factors A, and the matrix elements are the matrix 
elements taken in the Fermi sea multiplied by the amplitudes A. For instance. 



FS 



{ed n \A\ed U) - KX^yni^d Ti \A\ed U) 



FS 



{ed n \A\e ]di)^ AeATAitAt;(erf Ti \A\ed TDfs (A4) 

The energy can be obtained by summing the matrix elements of the Hamiltonian over all 
256 configurations divided by the norm of the states. However, since the central pairs (2 
and 3) and the shell pairs (1 and 2 , 3 and 4) are nor treated in a symmetric way, the cluster 
is not homogeneous. In order to force homogeneity, different amplitudes X^^^s and X^^^c can 
be defined for the shell and central pairs. Then equation (A4) has to be modified according 
to 

{ed n \ed Ti) - XlX,^X^^^Aed Ti \ed n)FS 
{ed Ti \A\ed Ti) - AL,,A^^X,.(e^ ti \A\ed n)FS 
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{ed n \A\e Ul)^ X,,,sXe^,Al^,X'^^^Aed U \A\ed n)Fs (A5) 

and in the minimization procedure we impose the constraint that the densities and correla- 
tions are the same everywhere in the cluster. Details can be found in Ref.0. 

Once we have the energy, the inverse compressibility is easily calculated as 

= ^ (A6) 

We define the point of phase separation to be the point where vanishes, i.e. where the 
system becomes unstable against local density fluctuations. A different definition of the 
phase separation boundary can be given by Maxwell construction. 

APPENDIX B: EFFECTIVE HAMILTONIAN FOR LARGE V 

In this appendix, we derive an effective Hamiltonian in the limit where the intersite 
repulsion V is much larger than the on-site interaction U and the hopping t. This can be 
done by using the Schrieffer- Wolff transformation^ 

= H + z[S,H]-^[S,[S,H]] + ... (Bl) 

We first split the Hamiltonian into two parts, H = H' + //mix) where H' mixes states within 
a given " Hubbard Band" , while -ffmix has matrix elements between different Hubbard bands 
typically separated by an energy of order U from each other. We can eliminate Hmix by 
choosing S so that 



H^,, = -t[S,H'] (B2) 



Then to first order in S 



H,s = H' + t^[S,H,,,,]. (B3) 

Using standard procedures, the resulting effective Hamiltonian we obtained consists of 
six parts 

20 



HeS — Hq + Hi + H2 + H3 + H4 + 



(B4) 



The leading term Hq — VH'V reads 



(B5) 



where 



(B6) 



i,a,<T' 



This is a Hubbard model with the additional constraint that two fermions cannot be on 
neighbouring sites. 

The remaining H^s are of order —t^jV. Hi describes the hopping of pairs 



Hi = -2^-Y,dld,+hx 



(B7) 



where dj = cj.T^,; i^i = Ci^^Ci^^ ) is the pair creation (annihilation) operator. Annihilation 
(creation) of pairs into (from) opposite spin fermions two sites apart is described by 



H2--yY. (4+i,T4-i,i + 4-i,T4+i,i) (^i+i + '^di + + h.c. 



H^ describes a pair-electron interchange 



(B8) 



while Hi stands for a simultaneous hopping of two fermions 



Finally, we have a diagonal term 



H.^ 



(B9) 



(BIO) 



(Bll) 



21 



REFERENCES 

* On leave from: Research Institute for Solid State Physics, Budapest, Hungary. 

Solyom, Adv. Phys. 28, 201 (1979). 

^ V. Emery, in Highly Conducting One-Dimensional Solids, edited by J. T. Devrese et al 
(Plenum, New York, 1979), p. 327. 

3E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 (1968). 

''p. Wzietek, F. Creuzet, C. Bourbonnais, D. Jerome, K. Bechgaard and P. Batail, J. 
Physique I (France), in press. 

^B. Dardel et al, Phys. Rev. Lett. 67, 3144 (1991); B. Dardel, D. Malterre, M. Grioni, P. 
Weibel, Y. Baer, J. Voit, D. Jerome, preprint. 

^H. J. Schulz, Phys. Rev. Lett. 64, 2831 (1990); Int. J. Mod. Phys. B 5, 57 (1991). 

^A. Parola and S. Sorella, Phys. Rev. Lett. 64, 1831 (1990). 

8 A. A. Ovchinnikov, Sov. Phys. JETP 37, 176 (1973). 

^F. Mila and X. Zotos, Europhys. Lett. 24, 133 (1993) . 
^°B. Sutherland, Phys. Rev. B 12, 3795 (1975). 
lip. Schlottman, Phys. Rev. B 36, 5177 (1987). 

12 P. A. Bares and G. Blatter, Phys. Rev. Lett. 64, 2567 (1990). 

13 F. D. M. Haldane, Phys. Rev. Lett. 45, 1358 (1980); ibid. 47, 1840 (1981); J. Phys. C 14, 
2585 (1981). 

i^G. Spronken, R. JuUien and M. Avignon, Phys. Rev. B 24, 5356 (1981). 

1^ For the Conformal Field Theory point of view, see e.g. H. Frahm and V. E. Korepin, Phys. 
Rev. B 42, 10553 (1990) and references therein. 

22 



C. N. Yang and C. P. Yang, Phys. Rev. 150, 321, 327 (1966); ibid. 151, 258 (1966). 
M. Fowler and M. W. Puga, Phys. Rev. B 18, 421 (1978). 
M. Ogata et al, Phys. Rev. Lett. 66, 2388 (1991). 
W. Kohn, Phys. Rev. 133, A171 (1964). 

X. Zotos, P. Prelovsek and I. Sega, Phys. Rev. B 42, 8445 (1990). 
B. S. Shastry and B. Sutherland, Phys. Rev. Lett. 65, 243 (1990). 
V. Emery, Phys. Rev. Lett. 37, 107 (1976). 

P. A. Lee, T. M. Rice and R. A. Klemm, Phys. Rev. B 15, 2984 (1977). 

With the noticeable exception of the U — +oo model at quarter filling, where the transition 
is effectively a 2kp CDW instability of spinless fermions: V. J. Emery and C. Noguera, 
Phys. Rev. Lett. 60, 631 (1988). 

see e.g. H. Yokoyama and H. Shiba, J. Phys. Soc. Japan 59, 3669 (1991). 

An exception to that statement can be found in W. Metzner and D. VoUhardt, Phys. Rev 
B37, 7382 (1988). 

T. A. Kaplan, P. Horsch and P. Fuldc, Phys. Rev. Lett. 49, 889 (1982). 

P. Fazekas, Int. J. Modern Phys. B3, 1765 (1989), P. Fazekas, in Condensed Matter 
Theories, Vol. 6, Edited by S. Fantoni and S. Rosatti, Plenum Press, New York, 1991, p. 
279 

M. C. Gutzwiller, Phys. Rev. 137, A1726 (1965) 

J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 (1966) 



23 



FIGURES 



FIG. 1. Phase diagram in the weak couphng hmit. 

FIG. 2. Phase diagram of the U — V model at quarter-fiUing. The nature of the various phases 
is discussed in the main text. 

FIG. 3. Typical configuration for N particles on L sites in the V = +oo limit: a) N = L/2; 
b) iV = L/2 + 1; c) and d) = L/2 - 1. Dots stand for empty sites. 

FIG. 4. Ground state energy as a function of flux for two cases: a) U = —1 and V = 8 and b) 
U = and V = 12. The boundary conditions are periodic for L = 12 and antiperiodic for L = 8 
and L = 16. The curves having a minimum at ^ = are nearly indistinguishble for L = 8, 12 and 
16, and we have plotted a single set of points. 
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